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Abstract 

We first reobtain in a simpler way the Haldane fractional statistics at thermal 
equilibrium using an interpolation argument. We then show that the mean 
occupation number for fractional statistics is invariant to a group of duality 
transformations, a nonabelian subgroup of fractional linear group. 



1 Introduction 

In the two-dimensional systems the particles are not constrained to obey the 
Fermi or Bose statistics but they may have fractional statistics; two well known 
exemples are offered by the quasi-particles in the fractional quantum Hall effect 
and by anyons. Prompted by this, Haldane ^ stated a generalized exclusion 
principle that allows for intermediate statistics, though his derivation of the 
particle occupation number does not involve the spatial dimensionality and so, 
it seems to be valid for an arbitrary number of dimensions. Wu 2 discussed 
the thermodynamics of the fractional statistics particles (geons) in the ideal gas 
approximation. Further developments on the thermodynamical aspects of geons 
can be found in [3], 0]. In an akin paper jQ to the present one, it was shown that 
the distribution funtion for geons enjoys a duality property. Quite recently, a 
number of authors ^ , starting from an appropiate definition for indistinguisha- 
bility in the phase space, have obtained a nontrivial classical counterpart of 
Haldane's exclusion principle. 

The purpose of this Letter is first, to rederive, using a very simple inter- 
polating argument, the occupation number distribution for geons at thermal 
equilibrium. Then, based on this derivation, we show that the relation for the 
occupation number remains invariant under a certain number of dual transfor- 
mations that organize in a group and we discuss the meaning of these dualities. 
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2 Model 



A very long time ago, when the black body radiation formula was still under 
scrutiny, a simple argument based on interpolation has been put forward in [7] . 
The argument goes as follows: Let be e the mean energy for a radiation field, 
averaged over time in a fixed unit volume. Using the equality of time averages 
and phase spase averages one can write for the mean squares values Ae^ of the 
fluctuating energy that 

Empirically, the mean square value gets contributions from low frequencies 
(the Rayleigh- Jeans part of the spectrum), where the mean energy e — —kT; 
therefore 

kT'^=e' (2) 

It also gets contributions from high frequencies (the Wien part of the spectrum) , 
where the mean energy e — hi' exp (— |^); hence 

kT'^ = h'^e (3) 

With the underlying assumption that the two contributions act indepen- 
dently and are of equal weight, relation (1) reads as 

kT'^^eh. + s' (4) 

It is more convenient to use the dimensionless variables n = -, — = a that 
lead to the differential equation: 



da 

The solution is 



^n + n'^ (5) 



n = — or e = p — - — ^ (6) 

exp (a — ao) — 1 / e — 

This is, the Planck law, here with a = merely an integration constant. 

kT 

This derivation emphasises the wave-particle character of bosons. 

How this reasoning can be extended for other statistics? Note that this type 
of interpolation works equally well for fermions. However, for fermions, as long 
as £ < /z, the long wavelength (high temperature) is negative: 
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and accordingly, relation (5) must be modified to 



dn 9 
— —=n-n^ (8) 

da 

When we take into account that n < 1, we obtain, as solution, the Fermi 
distribution: 

n = , ^ , , (9) 

exp [a — ao) + 1 

Introducing a genuine interpolation parameter g, i.e. a parameter ranging in 
the [0, 1] interval, one can write in a single stroke, for both fermions and bosons: 

dn , ^ n , ^ 

-^=n+{l-2g)n^ (10) 

with the solution 

n = ; \ 7Z (H) 

exp (a - ao) - (1 - 25) ^ ' 

Note that this represents also, as a particular case, for 5 = ^, the Maxwell 
distribution function. The above interpolating formula has also been obtained 
in Refs.8. 

At this moment we may ask what other modifications must be brought to (9) 
if one would allow for intermediate statistics betweeen Bose and Fermi ones on 
one hand, relation (9) suggests a series expension in powers of mean occupation 
number n. This is what we might expect because an interpolating statistics 
presumely contains more parameters and its complete characterization needs 
higher moments and/or higher powers of the first moment. On the other hand, 
higher order terms in (9) must be multiplied, for sure, by the factor g (g — 1), 
that excludes for standard fermionic {g = 1) and bosonic {g = 0) statistics any 
extra terms which, obviously, do not occur in (10). Limiting ourselves to the 
next higher order term n'"* and using again the independence and equal weights 
assumption (with the exception of g {g — 1) factor), we arrive at 

-^=n+{l-2g)n^+g{g-l)n'' (12) 
da 

The differential equation can be solved by an elementary integration. To 
this purpose and for later reference we make the change of variable: 

w='-^ (13) 

n 

The differential equation transforms to 

dw 10(10 + 1) ^, 

(14) 



da w + g 
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The solution is 



w9 + 1)1-9 = e("-"«) (15) 

Formulas (13) and (15) together gives implicitely the mean occupation num- 
ber as function of inverse temperature. The physical meaning of the parameter 
g can be drawn from the zero temperature limit {a — > oo) : 

? (16) 

Therefore, ^ represents the maximum number of particles a quatum state 
can accommodate. Actually, formulas (13), (15) represent exactly the parti- 
cle distribution for geons at thermal equilibrium obtained by Haldane P from 
a counting argument and from equilibrium condition requirements while (16) 
expresses the generalized exclusion principle. 

We have unfolded the derivation starting from the mean square value of the 
energy just to remind an historical argument [7]. Of course we could have started 
from the mean square value of the mean number directly. It is an well-known 
fact that the fluctuations in the number of bosons is super-poissonian: 



An2 =n + rf (17) 
while for fermions is sub-poissonian: 



An2 =n-n^ (18) 

The advantage of the above derivation of fractional statistics over the stan- 
dard one consists in its simplicity. The shortcoming lies in the fact the formulas 
for n are valid only at equilibrium. 

What are the transformations that leave invariant the occupation number 
formulas invariant, i.e. formulas (13) and (15)? To be more specific, we look for 
transformations for w that, while preserving the substitution (13) intact, makes 
relation (15), or, under the same boundary condition, the differantial equation 
(14), form- invariant. An earlier derivation for such a tranformation [2| and a 
theoretical prejudice lead us to consider transformations of the form: 

^ aw + b 

w = — (19) 

cw + a 

They are known as transformations belonging to the modular group SL(2,G), 
or the fractional linear group, where G stands for C, R or Z depending on 
whether the coefficients a, 6, c, d are complex, real or integer, respectively. The 

group multiplication law may be represented by a 2x2 matrix ~ ^ ^ 

with unit determinant ad-cd=l. The theoretical prejudice alluded above refers 
to the fact the SL(2,Z) has been found relevant for the complex conductivities 
in the quantum Hall effect 
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We substitute Eq. (19) into the diflFerential equation (14) and look for the 
coefficients a, b, c and d that leave the equation in same form. Rather long 
but straightforward algebra shows that, excepting the identity, there are only 
five possible tranformations of type (19) that leave the geon distribution form- 
invariant. These transformations (tk) are: 





W 


9 


a 


h 


1 


1 


a 


w 


9 


9 




iw 


9 


a 


iw + i 


.9-1 


.9-1 


U 


w + 1 


5-1 


a 


w 


9 


9 


tA 


1 


1 


a 


w + 1 


1-.9 


1-9 




iw + i 
i 


1-5 


a 



The first transformation has also been remarked and derived in a different 
way in j^. The others are new. Note that in some transformations (t2 and 
is) the coefficients must be imaginary to comply with the requirement of unit 
determinant ad-ch=\. 

These transformations deserve the name of dualities because they extend the 
theory from weak values of the statistical parameter g G [0,1] to strong values 
of the parameter g € [l,oo] or the other way aroimd. Indeed, starting with a 
weak g S [0, 1], transformations ti and ti lead to formulas valid on the whole 
real axis for g and on the range 5 G [1, 00], respectively. Starting with a strong 
g £ [1,00], transformation ^3 leads to formulas valid on the range g € [0,1]. 
Transformations ^2 and map the domain of validity for g (both g and g must 
be necessarily positive) into themselves. The dual distributions hold true at a 
temperature scaled with a function of g, at energies below ( t2 and ^5, a > 0) 
or above Fermi level (^1,^2 and ^4, a < 0). 

A distribution with a greater than unit g supports two interpretations: 

(i) as a distribution for particles; a physical realization is the Hubbard model 
PI, where due to the infinite repulsion, two sites can be occupied by at most 
one electron. (g — 2): and 

(ii) as a distribution for vacancies. For instance, the generalized exclusion 
principle (16) states that ^ is the maximum number of particles a quantum state 

can accommodate; then i , where g is the dual of g under ti would represent 
the minimum number of states available for one particle. 

The relationships between a specific distribution and its dual are given below 

g-n{a,g)+g-n{a,g)^l; (20.1) 
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{g-l)-n (a, g) + {g - 1) ■ n {a, g) = 1; 



(20.2) 



{g-l)-n{a,g)+g-n{a,g) = l; (20.3) 
g-n{a,g) + {g-l)-n{a,g) = l; (20.4) 

n(5,5)+n(a,5) = 0; (20.5) 

The last relationship is somehow peculiar. Apparently, due to g — > 1 — g 
transformation, it seems to convert the distribution for fermions or near fermions 
(5 > ^) into the distribution for bosons or near bosons (g < ^) and vice-versa, 
but the dual n {a,g) of n {a,g) is, in this case, negative 

The above dualities impose constraints on the expansion coefficients of ther- 
modynamic functions. For instance, the duality (20a) leads to a vanishing zero 
temperature heat capacity (Refs. 5-8). A detailed derivationand analysis of 
these constraints are relegated to an extended paper. Here we are content to re- 
mark that all five transformations are needed for these transformations to close 
upon themselves, i.e. to form a group. We think the group property is the most 
remarkable feature and we give below the table for the group composition law: 
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where to denotes the identity. The left mTiltiplication corresponds to entries 
along column, the right multiplic-atioii to entries along rows. It is worthy to 
note that transformations tj , j = 3.4 are inverse of each other, a fact that can 
also be seen from rels. (20.3), (20.4) 

3 Conclusions 

To sum up, we have proved that the fractional statistics distribution is invariant 
to a finite, nonabelian subgroup of the fractional linear group. 

We have discussed the duality in terms of transformations that, while leaving 
the theory invariant, relate two different "coupling constant" g regimes. There 
exists, also, another meaning of duality in this letter, as illustrated, for instance, 
in relation (4), namely, the wave-particle duality. The first term in the RHS 
of rel. (4) features the particle-like character of bosons and the second, the 
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wave-like character. The fact the two terms simply add togather expresses the 
complementarity of the two features. In standard quantum mechanics, this 
duality is embodied in the usual momentum-position commutator. We may 
wonder if the newly introduced term g {g — 1) for fractional statistics, in 
relation (12) signals the existence, in an analogous way, of another level of 
duality. 

In view of the above remarks it is tempting to put forward the following 
conjecture: there must exist another type of duality, in the sense above, but be- 
tween vortices and particles (as in some models for anyons [11]). The conjecture 
is rather speculative at the moment. We hope to return with a more rigorous 
statement in the next future. 
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